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ABSTRACT 

Aims. We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical 
properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy 
closest to the Milky Way, and compare these features with H i catalogs from the ATCA+Parkes H i survey. 

Methods. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distri- 
bution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy dis- 
tribution fits to 250, 350, and 500 jim emission from SPIRE observations using dust grain size distributions for graphite/silicate in low-metallicity 
extragalactic environments. 

Results. The dust cloud mass spectrum follows a power law distribution with an exponent of y = -1.8 for clumps larger than 4xl0 2 M Q and is 
similar to the Hi mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation. 
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1. Introduction 

The overall structure and evolution of the ISM of a galaxy is de- 
termined by the physical and dynamical conditions in the diffuse 
ISM where most of the processing of gas and dust clouds occurs. 
Thus inspecting the cloud properties by studying the atomic gas 
and the dust emission under different environmental conditions 
will help us to understand how this processing of matter pro- 
ceeds. 

In the present paper, we study the dynamical conditions of 
the diffuse gas in the Large Magellanic Cloud (LMC). We ex- 
amine optically thin dust emission at far-infrared (FIR) and sub- 
millimeter wavelengths by studying the spectral energy distribu- 
tions (SEDs) inferred from radiative transfer models. By com- 
bining an Hi synthesis survey and an imaging survey in the 
FIR and submillimeter emission performed by Herschel Space 
Observatory (Pilbratt et al. 2010), we investigate the detailed re- 
lationship between the diffuse atomic gas and interstellar dust 
emission in the LMC. The Herschel Space Observatory is de- 
signed to study the submillimeter and FIR emission from our 
universe and offers unprecedented resolution. The Herschel has 
a pay load of two cameras, PACS (Poglitsch et al. 2010) and 
SPIRE (Griffin et al. 2010). PACS employs two bolometer ar- 
rays for imaging, and SPIRE comprises spider-web bolometers 
with NTD Ge temperature sensors. Mapping strategy and data 
processing procedures are described in detail by Meixner et al. 
(2010). 

The Herschel survey of the Magellanic Clouds will provide 
information about the interstellar dust emission (Meixner et al. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



2006) from all phases of the ISM in the Magellanic Clouds. 
For diffuse atomic gas, we use data from the H i aperture syn- 
thesis survey in Kim et al. (1998) and data from the combined 
survey using single-dish observations with the Parkes 64-m ra- 
dio telescope (Staveley-Smith et al. 2003) that was published in 
Kim et al. (2003). The SED models provide spectral informa- 
tion in the 160 jt/m to 500 //m wavelength range. Relatively op- 
tically thin dust emission at FIR and submillimeter wavelengths 
can also provide information about the gas mass in the ISM be- 
cause the gas and dust are well mixed in most of the ISM phases. 
Therefore, the dust emission at the FIR and submillimeter wave- 
lengths can trace the column density and the structure of the 
ISM. 



2. Cold dust associated with atomic hydrogen 

2.1. Dust clumps 

The unprecedented quality of the SPIRE images improves 
our knowledge of the spatial distribution of galactic in- 
frared/submillimeter dust radiation and characterizes it as cloud 
features with clumps and filaments. We adopt the 'clump' ter- 
minology to represent entities with properties of typical 'clouds' 
seen in the molecular clouds in the Galaxy (Bergin and Tafalla 
2007). We examine the brightness distribution of the pixels treat- 
ing them as clumps by searching for peaks of local emission 
and adding the neighboring pixels of one clump to the objects 
(Williams, de Geus, & Blitz 1994). We find 7449 dust clumps 
generated by an automatic clump-finding routine for the 500 fim 
emission. The threshold being set to three times of an RMS of 
~0.3 MJy/sr (Fig. 1). The automatic clump-finding algorithm 
determines structure by first contouring the data at a multiple 
of the RMS noise of the data, searching for peaks of emission 
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Fig. 1. Dust clouds seen in the SPIRE 500 /jm image were 
identified using an automatic finding algorithm in which the 
lowest threshold was set to be 3 cr (blue contours), and they are 
overlaid on the Hi aperture synthesis image (Kim et al. 2003) 
of the LMC. la- is set to 0.3 MJy/sr at 500 /mi. The black box 
indicates the observed region in the Science Demonstration 
Program of Herschel Space Observatory (Meixner et al. 2010). 
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Fig. 2. Dust clouds seen in the SPIRE 350 /mi image were 
identified using an automatic finding algorithm, the lowest 
threshold being set to 3 cr (red contours), and overlaid on the Ha 
image of the LMC (Kim et al. 1999). The black box indicates 
the observed region in the Science Demonstration Program of 
Herschel Space Observatory (Meixner et al. 2010). 



that are local maxima in the SPIRE images, and then connect- 
ing pixels at each contour level from the highest to lower in- 
tensities. Isolated contours at each level are identified as clumps 
(Williams, de Geus, & Blitz 1994). 

By applying the same method to the 350 /mi image of the 
LMC, we find 8460 dust clumps (Fig. 2). More faint clumps are 
seen at 350 /mi than 500 /mi above the noise level. However, the 
dust clumps at 500 /mi are well correlated with the H i clouds 
and filaments. In general, the dust clumps at 500 //m identified 



in the present study are distributed quite uniformly in the LMC. 
Thus, for the clump analysis pursued in this study we use the 
clumps identified in the 500 //m image. 

For a clump consisting of a set of pixels with positions {x, y} 
and intensity {/}, the size of each clump is calculated based on its 
extent in the spatial dimension below. The clumps range in size 
from 9.8 to 47+1 pc with a median of 15+1 pc in radius given 
by 
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2.2. Physical properties of dust clumps 

Each dust clump was characterized by the radiation source 
properties associated with the dust clump and their ambient 
medium using DUSTY radiative transfer calculations performed 
by Ivezic et al. (1999). This model assumes a spherical geom- 
etry and supports an analytical form of the dust density distri- 
bution (Sarkar et al. 2006). The physical parameters constrained 
by fitting the observed SEDs include the chemical composition 
of the dust, grain size distribution, dust temperature at the inner 
and outer boundary, dust density distribution, optical depth at a 
specific wavelength, and the ambient interstellar radiation field 
(ISRF). We adopted the chemical properties of the dust deter- 
mined by Draine & Li (2001) and Weingartner & Draine (2001). 
For dust in the neutral medium, we adopted the size distribu- 
tion of Kim, Martin, & Hendry (1994, KMH) using the power 
law distribution function for grain sizes, n(d) oc cr q e~ a l cla with 
<7=3.5, a, m „=0.005 /mi, and a mflv =0.25 /mi. The classical model 
of grain size distribution was constructed by Mathis, Rumpl, & 
Nordsieck (1977, MRN), and Draine and Lee (1984) revised the 
MRN model with dielectric functions for graphite and silicate 
(Draine and Li 2001). 

Interstellar grain temperatures are calculated from 
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where B A {T) is the Planck function for dust at temperature T, 
Fj is the flux of the interstellar radiation field, and Qm{a, A) and 
Quv(fl,A) indicate the infrared and ultraviolet absorption effi- 
ciencies, respectively, and are functions of the dust grain radius 
a (Spitzer 1978). The ISRF is produced by three components: 
(1) emission from stars, (2) dust emission in the far-IR associ- 
ated with the heating of interstellar matter by absorption of star 
light, and (3) the cosmic background radiation (Cox et al. 1986; 
Chi & Wolfendale 1991; Strong et al. 2000). Unlike the Milky 
Way, confusion along the line-of-sight to the LMC is negligible 
because of its low inclination (Weinberg and Nikolaev 2001). 
The LMC also has low internal reddening, and its proximity per- 
mits the detailed study of its stellar population. Therefore, the 
ISRF in the LMC can be estimated with confidence. 

The dust equilibrium temperatures are obtained by deter- 
mining the energy balance between absorption and emission. In 
Fig. 3, we present three examples of applying DUSTY to dif- 
ferent clumps. The result of each fitting provides a characteris- 
tic dust emission temperature, Tj for each clump and an esti- 
mate of the graphite and silicate abundance ratio. The resulting 
graphite and silicate abundances from fitting the observed SEDs 
differ by 50% among the clumps in different environments in the 
LMC. The second approach to characterize clump temperature 
is to apply GRASIL (Silva et al. 1998). This spectrophotometric 
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self-consistent model computes the absorption and emission by 
dust in three different environments such as molecular clouds, 
diffuse ISM, and AGB envelopes with a dust model consisting 
of big grains, small grains, and polycyclic aromatic hydrocar- 
bons (PAH) molecules. The consistency check on the equilib- 
rium temperature is made using the GRASIL model and the ob- 
served SEDs (Fig. 3). In general, the observed fluxes at 160 - 
500 fim are ascribed to the dust in the diffuse medium heated by 
the interstellar radiation field. Heating of diffuse gas in the ISM 
can occur by means of the photoelectric heating in interstellar 
clouds where the dust grains act as catalytic surfaces. The resul- 
tant temperature from the fit to the observed SEDs ranges from 
15+0.4 to 25 ±0.4 K and is similar to the SPIRE temperature map 
(Gordon et al. 2010). 

At a given grain temperature, the mass of each dust clump is 
calculated from M d = S a><D 2 / K^xB^Td), where S a is the flux 
density at 500 //m and D is the distance of the galaxy at 50 kpc 
(Schaefer 2008), and B^iTd) is the value of the Planck function at 
500 /im, and a function of T d . We adopt a mean temperature for 
each clump as a function of the three categories shown in Fig. 3. 
The YSOs region and bar region in Fig. 3 indicate dynamically 
hot regions described in Kim et al. (1998, 2007), where k,\ is the 
mass absorption coefficient at 500 /im from the SPIRE observa- 
tions and a-,1 ~ 1.15 cm 2 /g is used. The derived dust masses are 
in the range of 1.8 x 10 1 M Q < M dust < 7.9 xlO 3 M . 

To characterize the clump properties, the mass spectrum of 
the clumps is examined. We choose the cumulative spectrum 
rather than the differential mass distribution since the relative 
masses will be more or less the same, when the mass range is 
small over the entire spectrum. For each catalog, we fit a power 
law function to the cumulative mass distribution (Rosolowsky et 
al. 2005; Kim et al. 2007): 



M V +l 



(3) 



where A^o is the number of clouds in the derived distribution with 
masses higher than 10 M and y is the index of the differen- 
tial mass distribution. The cumulative mass distribution of dust 
cores in the catalog of sources identified in this study and the 
residuals of the fit are presented in Fig. 4. The clump mass spec- 
trum follows a broken power law in the mass distribution. For 
masses lower than 4 x 10 2 M Q , it follows a power law distri- 
bution and decreases slowly with an exponent of the mass dis- 
tribution, y - -0.8. For the clumps more massive than 4 xlO 2 
M Q , the mass spectrum is significantly steeper with a slope of 
y - -1.8 + 0.1 than in the lower mass regime. The observed 
clump mass spectrum for 1.8 x 10 1 M < M dus , < 4.0 xlO 2 M , 
exhibits a flatter power law than the stellar initial mass function 
(Salpeter 1955). This confirms that most of the clump mass re- 
sides in massive clumps, while the most of the stellar mass is 
in low-mass stars (McKee and Williams 1997; Johnstone et al. 
2000). 

2.3. Gas and dust mass spectrum 

Observation and analysis of the dust clumps and neutral cloud 
structure indicate that a number of smaller clumps reside in the 
neutral diffuse clouds. The physical properties of dust clumps 
were characterized in terms of their mass spectrum, which fol- 
lows a power law distribution. We note that the distribution of 
dust clump mass is quite similar to the H i mass distribution, the 
index of the mass distribution being, y = -1.88. The mass of 
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Fig. 3. DUSTY model spectra and the best SED fitting results 
for three representative dust clumps in different environments 
(Top). The ISRF for each clump is 1.0 G (ISM region), 3.5 G 
(Bar region), and 7.0 G (YSOs region). GRASIL (Silva et al. 
1998) model spectra and the observed SEDs (Bottom). The data 
points are from IRAC, MIPS, and SPIRE. 



the H i cloud was determined from its integrated intensity / mea- 
sured in Kkm s 



M doud = 1.822 x 10 I8 /cr d m H , 



(4) 



where <x d is the total cloud area in cm 2 . Since the dust emis- 
sion in the far-IR continuum emission provides an estimate of 
the star formation rate (SFR) and grains remain strongly cou- 
pled to the neutral gas, we expect the mass spectra of the dust 
clouds to behave in a similar way to the mass spectra resulting 
in the diffuse gas. 

The mass function near the high-end appears to be even 
steeper with y - -2.1. The index of that power law fit can be 
steeper where star formation is ongoing, such that cloud dissipa- 
tion occurs (Wada et al. 2000). The similarity of clump mass dis- 
tribution suggests that the overall gas-to-dust mass ratio is more 
or less uniform in the region considered here. This is consistent 
with the overall gas-to-dust mass ratio being relatively uniform 
across the LMC, except for the super-shell affected by the su- 
pernova shocks (Gordon et al. 2003). A detailed analysis of the 
gas-to-dust mass ratio of two molecular clouds in the LMC can 
be found in Roman-Duval et al. (2010). 
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Fig. 4. Dust clump mass spectrum (Top). The slope of the mass 
spectrum, y, is -1.8 with an error estimate of +0.1. Hi clump 
mass spectrum (Bottom). The slope of the H i mass spectrum, y, 
is -1.88 with an error estimate of ±0.1 for the fit in the selected 
mass range. 



3. Summary 

Dust clumps have been identified and cataloged in the Herschel 
SPIRE survey of the LMC using an automated cloud-finding al- 
gorithm. The distribution of cold dust clumps is remarkably sim- 
ilar to the H i clump mass distribution, sharing an index of mass 
distribution, y - -1.8. However, the dust clump mass spectrum 
in the lower mass regime follows a flatter power law than the 
Salpeter stellar IMF. 
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